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Abstract: 

The cross-section for ZZ — > ZZ with arbitrarily polarized Z bosons is calculated within 
the electroweak Standard Model including the complete 0{a) corrections. We show the 
numerical importance of the radiative corrections and elaborate its characteristic features. 
The treatment of the Higgs-boson resonance is discussed in different schemes including 
the S'-matrix-motivated pole scheme and the background-field method. The numerical 
accuracy of the equivalence theorem is investigated by comparing the cross-sections for 
purely longitudinal Z bosons obtained from the equivalence theorem and from the com- 
plete calculation. In this context the full 0{a) corrections are also confronted with the 
enhanced corrections of O (aM^/s^M^), which were frequently used in the literature. 
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1 Introduction 



Gauge-boson scattering provides a window into the heart of spontaneously broken 
gauge theories: the gauge-boson self-interactions and the scalar sector, which drives spon- 
taneous symmetry breaking. Therefore, such processes found continuous interest in the 
literature [ |I|, |^, ^, ^, |^, ^ 0, |] since the very first years of spontaneously broken gauge 
theories. Since lowest-order predictions for all gauge-boson scattering amplitudes involve 
only interactions between gauge and scalar bosons, the corresponding cross-sections de- 
pend very sensitively on the non-abelian and scalar sector of the underlying theory. This 
sensitivity is even enhanced for high-energetic, longitudinally polarized massive gauge 
bosons, owing to the presence of gauge cancellations. A longitudinal polarization vector 
contains a factor fc°/M, where k'^ and M are the momentum and mass of the correspond- 
ing gauge boson, respectively, and induces contributions to the matrix element that grow 
with energy. In spontaneously broken gauge theories such contributions cancel in the high- 
energy limit, as required by unitarity. For 't Hooft gauge- fixing conditions these so-called 
unitarity cancellations are quantitatively expressed by the Goldstone-boson equivalence 



theorem (ET) [ |I0| [Til, |T2|, ^ which relates amplitudes for longitudinal gauge bosons 
to those of the corresponding Goldstone bosons and thus reflects the connection between 
gauge and scalar sector of the theory. 

In the minimal SU(2) x U(l) electroweak Standard Model (SM) only one physical 
scalar field remains after spontaneous symmetry breaking, viz. the Higgs boson, which 
plays a central role in the discussion of massive gauge-boson scattering. Virtual Higgs- 
boson exchange is needed to prevent the 2 — > 2 scattering amplitudes of longitudinal 
gauge bosons from violating the (perturbative) unitarity bound at high energies. In 
turn, the requirement of unitarity can be used to derive bounds on the Higgs-boson mass 
below which the SM remains weakly interacting and treatable in low-order perturbation 
theory [ 0, 0. These bounds are of the order of 1 TeV and are slightly strengthened by 
including the O {aM^/ s'^M^) radiative corrections (RCs) to gauge-boson scattering [ 
H . As already pointed out in Ref . [ |^ , all these bounds are only qualitative, since they 
are obtained by applying perturbation theory in a region where it breaks down. The 
bounds on Mh can be related to a scale of new physics, which is necessary to avoid the 
Landau pole in the scalar self- interaction [|^. If these bounds are not satisfied the Higgs 
sector becomes strongly interacting. In this case large effects of new physics should arise 
in gauge-boson scattering and these processes would be particularly suited to study the 
electroweak symmetry breaking sector of the SM [ ^ |IU|, |H| . 

Gauge-boson scattering reactions can be studied at all high-energy colliders, i.e. pp col- 
liders like the LHC, e"'"e~ colliders like the NLC, or fi~^fi~ colliders, where these reactions 
naturally appear as subprocesses. At high energies {E 3> Mw) the incoming particles 
radiate plenty of gauge bosons. Similar to the well-known Weizsacker- Williams approxi- 
mation for photonic reactions also massive vector-boson scattering at high energies can be 
approximated by convoluting the vector-boson cross-section with the corresponding flux 
of gauge bosons. This approximation is known as equivalent vector-boson method (see e.g. 



Ref. [ |T^] and references therein) . 



At high energies, where the investigation of gauge-boson scattering is most interesting, 
the RCs are typically large and need to be taken into account. In this paper we investigate 
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the effects of RCs on on-shell massive gauge-boson scattering processes. We have chosen 
the simplest representative, the process ZZ — >■ ZZ. It contains all interesting features that 
are typical for massive gauge-boson scattering such as the occurrence of a Higgs-boson 
resonance or enhanced RCs associated with a heavy Higgs boson. On the other hand, 
complications by bremsstrahlung corrections, which occur for W bosons, are absent. 

We calculate the complete 0{a) RCs to ZZ ZZ and present a detailed numerical 
discussion of the C(a)-corrected cross-sections both for the unpolarized case and the most 
interesting individual polarizations. Once RCs are taken into account, the introduction 
of a finite decay width of the Higgs boson, which is necessary for a sensible description 
of the resonance, is non-trivial owing to problems with gauge invariance. We compare 
different treatments such as the naive introduction of a finite width, Laurent expansions 
about the complex pole, as well as Dyson summation of self-energy corrections. The latter 
procedure is, in particular, applied within the framework of the background-field method 
(BFM) (see Ref. [ |T6[ and references therein), where Dyson summation does not disturb 
the underlying Ward identities [ |l3l which guarantee gauge cancellations and unitarity. 

For longitudinal gauge-boson scattering the radiative corrections of O (aM^/s^M^), 
which dominate for a heavy Higgs boson, have been calculated in the literature using the 
ET. We test the accuracy of such an approach by comparing these results with the full 
0{a) corrections. Moreover, we have calculated the 0{a) corrections as predicted via the 
ET, which possess a very simple analytical form. 

This paper is organized as follows: After some preliminary remarks in Section ^ 
about kinematics, conventions, and discrete symmetries, we discuss the lowest-order cross- 
sections in Section ^ In Section ^ we describe the explicit calculation and the structure of 
the 0{a) corrections. The different methods for introducing a finite Higgs-boson width are 
presented in Section ||. A brief description of the apphcation of the ET to ZlZl — > ZlZl 
and the heavy-Higgs-boson effects in Section ^ concludes our presentation of the calcu- 
lational framework. Numerical results are discussed in Section |^, and Section ^ contains 
our conclusions. Appendix |^ provides a further discussion of the Landau singularity that 
appears in some box diagrams. In Appendix |^ we present the full analytical results for 
the 0{a) corrections obtained via the ET. 

2 Preliminaries 

2.1 Kinematics and conventions 

We consider the reaction 

Z{h, Ai) + Zih, X2) ^ Z{h, A3) + Z{k,, A4) , (1) 

where ki and Aj denote the momenta and helicities of the incoming and outgoing Z bosons, 
respectively. We use the indices L, T, and U to indicate longitudinal (A = 0), transverse 
(A = ±), and unpolarized Z bosons, respectively, and characterize definite polarization 
combinations by a sequence of four letters, e.g. LTLT stands for ZlZt —>■ ZlZt. 
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The incoming particles travel along the z axis and are scattered into the x-z plane. 
In the center-of-mass system (CMS) the momenta and polarization vectors Ei[\i) read 



K = 


(E, 0, 0, -V) , 




4(0) = 


(-P, 0, 0, E) /M^, 


4*(o) 


^^(±) = 


(0, -1, ±i, 0)7^2, 


4*(±) 




(E, 0, 0, p) , 




4(0) = 


(-P, 0, 0, -E)/Mz, 


4*(o) 


4(±) = 


(0, 1, ±i, 0)/v^, 


4*(±) 



(-E, — psin^^, 0, — pcosi 



(-E, psin6', 0, pcosO) , 
4'*(0) = (p, Esin^, 0, Ecos^) /Mz, 

(0, cos^, Ti, -sine)/y2 (2) 



in terms of the energy E of the Z bosons, their momentum p = \J E"^ — M|, and the 
scattering angle 9. The Mandelstam variables are defined as 

s = {ki + hy = AE\ 

t = {ki-ksf = -Ap^sm^e/2, 

u= {ki- kif = -V cos2 e/2 . (3) 



Following the treatment of Ref. [ [T^] for 77 WW, we introduce the 83 standard 
matrix elements (SMEs) Aiijki which contain the complete information about the boson 
polarizations^. The invariant matrix element Ai is decomposed into a linear combina- 
tion of the SMEs with invariant functions Eijki{s,t) as coefficients. Exploiting discrete 
symmetries, the number of independent SMEs can be reduced. 

In terms of the invariant matrix element Ai the differential cross-section is expressed 

as 

da\ 1 , . . ,2 



■. / A1A2A3A4 



l-^A,A.A3A4l • (4) 



The unpolarized cross-section results from an average over the initial states and a sum 
over the final states, 

\ / unpol Ai,A2 A3,A4 \ / A1A2A3A4 

More generally, the correct average is obtained by multiplying with 1/3 for each unpolar- 
ized Z boson and by 1/2 for each transverse Z boson in the initial state. 
The integrated cross-section is obtained by 

1 ..360° /•i80°-6»cut da 

c^tot = o / / d^ sin^ — , (6) 

2 Jo° J Scut dSz 

where 6'cut denotes an angular cut which is set to 10° in our numerical evaluations. The 
symmetry factor 1/2 results from the presence of two identical particles in the final state. 

^Only 81 SMEs are linearly independent, the other two are kept for convenience. 
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2.2 Discrete symmetries 

As a consequence of Bose symmetry the amplitude Ai is invariant under tlie inter- 
change {ki,ei) ^ {k2,e2) and/or {h,el) ^ (k^^el), i.e. 

Mx,x2x^iX4{E,6) = Mx2M\A\:i{E,6) 

= Mx,x,x,X:AE, 180° + 6)= MxM^xAE, 180° + 6). (7) 

This implies for the cross-sections 

\ / A1A2A3A4 \ / A2A1A4A3 

\ / A2A1A3A4 \ / A1A2A4A3 

In particular, all cross-sections with equally polarized incoming and/or outgoing Z bosons 
are forward-backward symmetric. 
CPT symmetry entails 

'^AiA2A3A4 = -^A3A4AiA2 (9) 

and the analogous relation for the cross-sections. 

Because quark mixing is completely negligible for this process, we use a unit quark- 
mixing matrix, and thus also CP is an exact symmetry^. As a consequence, the helicity 
amplitudes are related as follows 

A<AiA2A3A4 = -A^-Ai-A2-A3-A4, (10) 

and the cross-sections do not change if all helicities are reversed. 

Owing to CP invariance all SMEs involving the totally antisymmetric Levi-Civita 
tensor drop out, and only 43 SMEs can appear. As a consequence of Bose and CPT 
symmetry only the sum of each SME and the ones obtained from the interchanges 

ki, £3, fcs) ^ (£2, ^2, £4, ki) and (ei, fci, £2, ^2) ^ (£^3, ^3, ^X, ^4) occur. This leaves 
19 SMEs, among which 17 are independent. 

3 Lowest-order cross-section 

To lowest order, only the three diagrams in Fig. ^ contribute and yield the following 
amplitude: 

cl,sl,\s-Ml^ t-Ml^ u-Ml) ' ^''> 



where cw = cos^w = My^/Mz, sw = sin^w = y^l — c^, and the relevant SMEs J^qq are 
given by 

= (ei ■ e^M ■ el) > -^S = (^1 ■ ^;)(^2 ■ el) , = {e, ■ el){e, ■ el) . (12) 

Explicit formulae for TVlBorn for the 81 polarization combinations are listed in Table |l|. 

^Even for a general quark-mixing matrix, CP would be violated in the considered process only beyond 
the one-loop level. 
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Figure 1: Lowest-order diagrams for ZZ TIL 
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Table 1: Polarized lowest-order matrix elements {S — s — M^, T — t — M^, U — u — 
and s = s - 4M|) 
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Polarization 



27r 



X (entry) 



LLLL 



— cos e^cut 



16M| 



LLTT 



cos 6*, 



cut 



LTLT 



cos6', 



cut 



LLLT 



{Ml + 2M|)^ f 

2M|s \ 2M^ + s(l 



cos^cut) _ 2 cos gcut (6Mg + s sin^ gcut^ 



cos 6'cut^ 



AMI + s sin^ ^, 



cut 



LTTT 



^ <^21n 



s (1 + cos 



2M2 -COS^cut) 



cos^cut(20M2 +3ssin2 
4M2 + s sin^ e,^t 



TTTT 



-8 In 



s (1 + cos 6*, 



cut; 



2M2 +s(l-COS^eut) 



+ 



cos6'cuts(8 + llsin^ 6'cut) 



4M^ + s sin^ 6^, 



cut 



Table 2: Polarized lowest-order cross-sections at high energies (s ^ M|,M^) integrated 
over ^eut < ^ < 180° - ^eut 



The dimensionful ZZH coupling leads to a suppression of the Born matrix element by 
a factor M|/|s — for \s — ^ M|. As a consequence, the lowest-order matrix ele- 
ment for purely transverse Z bosons is suppressed for high energies, s ^ M^, by a factor 
M^/s. Each longitudinal Z boson introduces a factor y/s/Mz via its polarization vector. 
For helicity amplitudes with more than two external longitudinal Z bosons unitarity can- 
cellations take place such that the Born matrix element with purely longitudinal Z bosons 
approaches a constant and those with three longitudinal Z bosons behave as M^j \fs at 
high energies. As a remnant of the unitarity cancellations the matrix elements involving 
four and three longitudinal Z bosons are enhanced by a factor M^/M| if M| <C -C s. 

The analytical results for the asymptotic behavior of the integrated cross-sections 
at high energies (s ^ M|,M^) are listed in Table ^ The cross-section for purely 
longitudinal Z bosons (LLLL) and the ones with two transverse and two longitudinal 
Z bosons (LLTT, LTLT) behave as 1/s for high energies. The cross-sections with one 
or three longitudinal gauge bosons (LTTT, LLLT) are proportional to and the 

cross-section for purely transverse Z bosons (TTTT) vanishes as at high energies for 
^cut > 0. When integrated over the full scattering angle the TTTT cross-section behaves 
as owing to the t- and u-channel poles. The cross-sections not shown in Table |^ are 
obtained from the symmetry relations (|) and (P), e.g. ctlllt = ctlltl = '^cttlll, where 
the factor 2 originates from the different spin average. 

Figure |^ illustrates the (exactly calculated) polarized lowest-order cross-sections in- 
tegrated over 6'cut < 9 < 180° — 6'cut for ^cut = 10° and a low Higgs-boson mass of 
Mh = 100 GeV. In Fig. ^ we show the same cross-sections for Mh = 700 GeV using a 
naive constant finite width to describe the Higgs-boson resonance, as discussed in Section ^ 
below. The enhancement of the LLLL and LLLT cross-sections caused by the factor 
M^/M| can be seen by comparing Figs. |^ and |^. The Higgs-boson resonance occurs only 
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Figure 2: Integrated lowest-order cross-sections for various polarizations at Mh = 100 GeV 
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Figure 3: Integrated lowest-order cross-sections for various polarizations at Mh = 700 GeV 
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for equally polarized Z bosons in the initial and final state, Ai = A2, A3 = A4. At lowest 
order, the LLLL cross-section dominates independently of the Higgs-boson mass. 



4 Radiative corrections 

4.1 Calculational framework 

We have performed the calculation of the radiative corrections (RCs) in 't Hooft- 
Feynman gauge both in the conventional formalism and in the background-field formalism, 
applying the on-shell renormalization scheme in both cases. We follow the conventions 
of Ref. [ |18[ for the conventional formalism and of Ref. [ |l6l for the background-field for- 



malism. In the conventional formalism the field renormalization is fixed such that no 
external wave-function renormalization is needed. In the renormalization scheme intro- 
duced in Ref. [ |l6l for the background-field method the field renormalization is determined 



by gauge invariance, and a non-trivial external wave-function renormalization is required. 



as explicitly described in Ref. [ |T3 



The Feynman graphs have been generated and drawn with FeynArts [ |TP|. Both in 
the conventional formalism and in the background-field method we have performed two 
independent calculations. 



One evaluation is based on the calculational method described in Ref. [ |T8| . With 
the help of Mathematica [ ^ the amplitudes are decomposed into SMEs and invariant 
functions, and the one-loop contributions to the invariant functions are expressed in terms 
of standard tensor integrals. The tensor integrals are reduced to the standard scalar one- 
loop integrals, as described in Ref. [21]. The scalar one-loop integrals are evaluated 



using the methods and general results of Ref. [ ^ . The last two steps are performed 
numerically using own Fortran routines. 

In the other calculation the algebra is performed with Mathematica and Form [ ^ 
and has been partially checked with FeynCalc [ The resulting symbolic amplitudes 
are automatically converted into a Fortran program. Instead of using SMEs, all scalar 
products of four-vectors are grouped together and calculated at run-time by inserting 
the explicit representations (0) for the polarization vectors. The tensor integrals are 



numerically reduced to scalar integrals, which are evaluated using the FF package [ |25|. 
The code thus obtained executes favorably fast and numerically stable. 

Because of the length of the results we do not list the analytical expressions but give 
only an inventory of the 0{a) RCs and discuss some important features. 

4.2 Inventory of 0{oi) corrections 

Both in the conventional formalism and in the background-field formalism about 550 
Feynman diagrams contribute to ZZ — > ZZ at one-loop order. The one-loop corrections 
can be classified into self-energy corrections, vertex corrections, box corrections, and wave- 
function-renormalization corrections. All of them can be divided into s-, t-, and u-channel 
contributions, which are related by simple transformations. In the following we list only 
the s-channel Feynman graphs for the conventional formalism in 't Hooft-Feynman gauge. 

The diagrams contributing to the self-energy corrections in the s channel are shown in 
Fig. H. The diagram with a virtual W and field actually represents two diagrams with 
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Figure 4: s-channel self-energy diagrams 




Figure 5: s-channel final-state vertex diagrams 

opposite charge flow. After renormalization the diagrams of Fig. ^ yield together with 
the corresponding t- and w-channel diagrams the following contribution to the invariant 
matrix element: 

where is the renormalized Higgs-boson self-energy. 

For each of the six vertices appearing in the Born diagrams (Fig. [l|) there is a set of 
vertex corrections. In Fig. ^ we show the diagrams that constitute the corrections to the 
final-state vertex in the s-channel diagram of Fig. 0. Note that each graph with three 
charged fields or two different charged fields in the loop represents two diagrams with 
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opposite orientation of the charge flow. Owing to the simple tensor structure of the ZZH 
vertex, the vertex corrections have the relatively simple form 



5-M.ertex = -^ E ^772 [2MS^o"^''W + (-^S + M;yr^(r)] (14) 

with the two renormalized form factors F^^^ir^ and F^'^^ir) for each channel, the 
corresponding SMEs A^qo^ from (0), and 

Mtl = (ei ■ e2){el ■ h){e: ■ h)/s, A^g = (e* ■ ^^(^i ■ k2){e2 ■ h)/s, 

Mi;i> = {ei ■ el){e2 ■ h){e; ■ k2)/s, M^J = (£2 ■ e*,){ei ■ k,){el ■ k,)/s. (15) 

The factors 1/s in (p!5| ) have been introduced to render the matrix elements dimensionless. 

The s-channel box diagrams (i.e. those with natural variables s and t) are shown in 
Fig. ^. Note that again all graphs with three or four charged fields in the loop represent 
two Feynman diagrams with opposite orientation of the charge flow. The analytical ex- 
pressions for the box diagrams are rather involved and require all CP-conserving SMEs. 
It turns out that the results for the bosonic box diagrams are shorter by a factor of 
about 3/2 in the background- field formalism as compared to the conventional formalism. 
Because of the involved structure of the fermion-Z-boson couplings, the most compli- 
cated expressions are those for the fermionic box diagrams, which are identical in both 
formalisms. 



Following the complete on-shell renormalization scheme of Ref. [ |T8| for the conven- 



tional formalism the field renormalization is chosen such that no extra wave-function 



renormalization is necessary. However, in the on-shell renormalization scheme of Ref. [ |T| 
for the background-field formalism a non- vanishing wave-function renormalization for ex- 
ternal particles is required. In a strict 0{a) calculation the wave-function renormalization 
corrections are given by 

5M„i = 2 5RzM Born; 

(16) 

where 

5Rz = -Re{j:'^^iMi)} (17) 

is the 0{a) contribution to the wave-function renormalization constant Rz = 1 + SRz- 
The function E^^(s) = dSf^(s)/ds denotes the derivative of the transverse part of the 
renormalized Z-boson self-energy. 

Thus, the full one-loop matrix element reads 

5Ali.loop = SMsclf + vertex + ^Aibox + SM^f. (18) 

4.3 Corrected cross-section 

It turns out that the 0{a) corrections are comparable or even larger than the lowest- 
order contributions for various important configurations. In order to obtain meaningful 
predictions, it is therefore necessary to consider not only the interference between the 
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Figure 6: s-channel box diagrams 
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lowest-order and the one-loop matrix element but to take into account the complete 
square of the matrix element and to define the corrected cross-section as 

^ = -KJ-r l-^Born + 57Wl-loop|' (19) 

including the terms involving |5A^i-ioop|^- 

In this way we end up with 0{a) accuracy where A^Bom dominates and 0{1) accuracy 
otherwise (relative to the leading loop order). To obtain 0{a) accuracy everywhere the 
(9(a^) corrections would be required. Note that the interference between 5M.2-\oop and 
-^Born is supprcsscd with respect to |(5A1i_ioopP if the importance of 5A^i-ioop results from 
a suppression of A^Bom and not from the presence of an effective large loop-expansion 
parameter. 

In Figs. 1^ and |^ we show the corrected polarized cross-sections as defined in (plQl) . 
The Higgs-boson resonance in Fig. ^ is treated by Dyson summation within the BFM [cf. 
Section p?2| , (p6D]. In the case of purely transverse Z bosons (TTTT) the cross-section is 
drastically enhanced compared to the lowest order and behaves as 1/s at high energies 
(without cut-off it even would not go down with energy). For a small Higgs-boson mass 
this cross-section becomes the dominating one. Apart from the LTTT case, the corrections 
reach the size of the lowest-order cross-sections for all polarizations. This is probably due 
to corrections of the form0 (a/7r)(ln(s/M|))^ 0.2 that are further enhanced by numerical 
factors and t- and w-channel poles. The relative corrections depend only weakly on the 
Higgs-boson mass apart from the polarizations LLLL and LLLT, where the corrections 
involve extra factors M^/Mz- As a consequence the LLLL cross-section dominates for a 
large Higgs-boson mass. As for the lowest-order cross-sections, the Higgs-boson resonance 
only contributes for Ai = A2 and A3 = A4. For the TTTT channel the resonance is 
proportional to the corresponding strongly suppressed Born cross-section and thus not 
visible. 

Because the lowest-order cross-sections are not dominating, the universal corrections 
associated with the running of a and the p parameter, which are related to the lowest 
order, are not leading. 

4.4 Landau singularities in four-point functions 

The four-point function [ Q 

Do(M|, M|, M|, M|, t, u, m^, m^, m^) , (20) 
exhibits a Landau singularity [ of the form 

= - (21) 
/A - ie 



for t < 0, u < 0, and Mz > 2m, where D^q^ is regular and 



tu 

A = — 
16 



t - 4m2) (u - 4m2) - (2M| - Am^ 



(22) 



■^At high energies vertex and box corrections typicaUy yield contributions of this kind, as e.g. exphcitly 
calculated for e+e" W+W" in Ref. [ |2§. 
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Figure 7: Integrated corrected cross-sections for various polarizations at Mh = 100 GeV 
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Figure 8: Integrated corrected cross-sections for various polarizations at Mh = 700 GeV 
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Figure 9: The location of the Landau singularities in the Dq function (^) for m = 0. The 
dashed lines at ^cut = 10° and 180° — 6'cut = 170° indicate the integration interval. 

With b = M| — m^s and = E'^ — M| this becomes 

A = {p^ sin^ e - b)p^ sin^ 6 . (23) 

Squaring the matrix element promotes the root singularity at p^ sin to a pole 

which is not integrable and thus leads to a formally divergent cross-section. 

This singularity should disappear from physical observables. The condition Mz > 2m 
suggests that it is related to the instability of the Z bosons. In fact, as illustrated in 
App. 0, it is canceled by diagrams that contribute to the inclusive process ZZ 4/, 
which cannot be separated from ZZ — >■ ZZ once the decay of the Z bosons is taken into 
account. Moreover, one should notice that colliding Z bosons which are radiated off from 
incoming particles possess an invariant mass < so that the condition > 4m^ is never 
fulfilled in the physical region of phase space. The use of on-shell Z bosons (g^ = M|) is 
just part of the equivalent vector-boson approximation. 

The Landau singularity appears in practice for box diagrams involving light fermions, 
i.e. with m = -C Mz- The location of the singularity in phase space for m = is 
shown in Fig. |. It appears at 9 = 90° for p^ = M|, i.e. = 2V2Mz ^ 258 GeV, and 
moves fast towards the forward and backward directions with increasing energy. Its effect 
is most prominent at low energies and becomes small at high energies. Moreover, it is 
located outside the angular region 10° < 9 < 170° for ^/s > 500 GeV for all fermions. 
In the following we always consider the cross-section in regions where this singularity is 
absent or negligible. 

5 Higgs-boson resonance 

Diagrams that involve a Higgs-boson propagator in the s channel have a pole at 
s = My. The double pole in the self-energy diagrams of Fig. ^ is reduced to a single pole 
after Dyson summation. 
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If Mh > 2Mz a proper treatment of the resonance is necessary to render the cross- 
section finite and meaningful. The naive introduction of a finite Higgs-boson width via 
the substitution 

(24) 



s - s - + iMnrn 



or naive Dyson summation 



(25) 



amounts to an inclusion of an incomplete set of higher-order corrections such that the 
resulting matrix element becomes gauge-dependent and violates the Ward identities and 
thus also the gauge cancellations, which guarantee unitarity. 

5.1 Pole expansion 

Since the poles of the S matrix are gauge- independent, it has been proposed [ ^ to 
perform a Laurent expansion about the complex pole. In a naive way this can be done by 
decomposing the contributions of the resonant diagrams into resonant and non-resonant 
parts and introducing the finite width only in the former. For vertex corrections this leads 
to the substitution 



vertex 



„2 2 

<-W^W i=0,l 
„2 „2 

^W-^W i=0,l 



is) 



+ M 



(s)^ 
Oi J 



7ZZH1 



(s)^ 
Oi , 



+ 



'ZZH 



{Ml) 



-Ml 



(26) 



where [to 0{a) accuracy] 



rH = ImE^(M2)/MH 



(27) 



is the decay width of the Higgs boson. For the lowest-order and self-energy contributions 
we write 



„2 2 



M 



00 



s-Ml 
1 



C\\tS 



'00 



WW 



Ml 



1 - 



s-Ml 

S^(Mg) 
s-Ml 



S'^(M^) 



s 

S^(.)-S^(M^)-(s-M^)S'^(M^) 



Ml 



e^Mj 



M 



00 



1 



(1 - s'^(M^; 



S-MI + iMnPH 
S^(.)-S^(Mg)-(.-Mg)S-^(Mg) 
is - Miy 



(2^ 



with Ti'^{s) = dS^(s)/ds. We have used the fact that the renormalized Higgs-boson 
self-energy fulfills ReS^(M^) = and that ImS^(M^) = MnTn is accounted for by the 
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resummed terms. As the inclusion of a finite width corresponds just to the inclusion of 
higher-order terms, all expressions are equivalent at the one-loop level. 

According to (^), Fjy is a tree-level quantity if is calculated at one- loop order. 
Consequently, in a one- loop calculation one ends up with accuracy on resonance and 
different treatments of the finite Z-boson width differ relatively in 0{a). An improvement 
to 0{a) accuracy on resonance would require to use the fully C(a)-corrected Higgs-boson 
width in the resonant denominators. 

The above procedure is gauge-independent because we modify the amplitude by terms 
that depend only on the gauge-independent residue of the pole and the physical mass 
and width. However, the actual application of the pole expansion deserves some care. 
Firstly, the above treatment is not uniquely determined by the resonance pole, because 
the Laurent expansion is only applied to the form factors but not to the SMEs and the 
split-up between these two is not unique. On the other hand, the terms introduced by the 
modification of the amplitude in general violate the Ward identities and thus eventually 
unitarity. This problem could be avoided and the pole scheme could be uniquely defined by 
including the complete matrix elements into the Laurent expansion. This, however, leads 
to problems in defining the residues, i.e. in particular the corresponding momenta and 
wave functions, for more general processes in certain kinematical regions (cf. Ref. [ 
In the following we show how one can exploit the above-mentioned freedom in the pole 
expansion in order to eliminate unitarity-violating terms. 

We first illustrate the procedure at tree level. A general pole expansion is obtained by 
absorbing some arbitrary function f{s) with /(M^) = 1 into the SME M.'qq and perform- 
ing the Laurent expansion for the resulting modified form factor. After resubstituting the 
original SME this amounts to the replacement 

(29) 

The added terms are proportional to f{s)M'o'f^Tn/{s - M'^)/{s - M'^ + iMnrn). If the 
matrix element f{s)A4'^Q grows too fast with energy these terms violate unitarity at high 
energies. In the high-energy limit the ratio between the Born cross-sections for longitu- 
dinal Z bosons (LLLL) calculated for Fh 7^ and Fh = behaves as 1 + /(s)^F^/9M^ 
for real f{s), i.e. f{s) = 1 yields a result that is off by a constant factor, but for instance 
f{s) = M^/s reproduces the correct high-energy limit. While different choices of f{s) 
by construction do not modify the resonant contribution they differ evidently in the non- 
resonant terms. This indefiniteness of the non-resonant lowest-order contributions gives 
rise to ambiguities of relative 0{a) in the resonance region. 

At the one-loop level the generalized pole expansion is obtained by absorbing arbitrary 
functions fij{s) with fij{M^) = 1 into the SMEs Aiif before performing the Laurent 
expansion of the form factors. Besides the appearance of several functions fij such a 
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2 2 ■'^'OO 



fis) 



S-MI + IMhFh 



1 - m 



general pole expansion even includes terms involving their derivatives f-y For our purposes 
it is sufficient to consider the following modified pole expansion 



32 Ml 



vertex 



„2 2 



E(Mo 

! = 0,1 



:^ZZH 



{Ml)f{s) 



:^ZZHl 



s-Mi + iMnrn 



+ 



Born 



self 



M 



s - + iMHlH 



s)-FP"{Ml)f{s) 
)S'^(M^)/(.) 



(30) 



which differs from a consistent expansion with = f{s) only by terms of the order 



of r^, i.e. O^a^). If we introduce the finite width as in p6| ) and (pH]), i.e. with /(s) = 1, 
we modify the cross-section for longitudinal gauge bosons at high energies by a constant 
contribution of 0{aT-ii/MY{). In the modified version of the pole scheme the high-energy 
behavior can be improved by choosing a suitable function /(s) that vanishes sufficiently 
fast at high energies. With our definition of the SMEs (|r^) it is sufficient to choose /(s) = 
M^/s. Note that if we did not include factors M^/s to render the SMEs dimensionless 
we would obtain a contribution to the matrix element that grows with s for f{s) = 1, i.e. 
that violates unitarity. The freedom parametrized by f{s) in (|30|) affects the non-resonant 
contributions in 0{a). On resonance this introduces ambiguities in Ola"^) relative to the 
leading resonant terms. 

The above recipe for the usual on-shell renormalization scheme is directly connected 
with an expansion of the transition matrix element about its complex pole. In Ref. [ 



such an expansion was explicitly described for angular-independent resonances, where the 
complications [|2^ in defining wave functions and momenta on resonance are absent. The 
procedure of Ref. [ ^ can be directly transfered to ZZ TIL. In this respect only the 
angular-independent, i.e. the (one-particle-reducible) lowest-order, self-energy, and vertex 
contributions in the s-channel, A^ipr, are relevant. After Dyson summation these can be 
written (assuming no truncation of the perturbation series) as follows: 



E 



WW jj=o,i 



(^o. + i^^^^(^))(^o.+i^-^^^(^)) 
s - M2 + S^(s) 



(31) 



product of full vertex functions and the full propagator. The additional SME 
M is defined as 



(£1 ■ k2){e2 ■ k\){el ■ k^){e\ ■ /cg)/^^ 



^11 

The complex pole Sp of ( |3TD is obtained as the solution of 



Sr, 



0. 



(32) 



(33) 



(s) 

Since Ai\p^ is analytical in s, it can be continued to complex s and expanded about s 
The leading term in this expansion is given by the resonant part 



Mr 



(s) 



(s-s,)(l + S'^(.,)) 



(34) 
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The residue of A^reso at the pole — Sp) can be interpreted as the product of two 
physical am] 
replaced by 



physical amplitudes M.^^^^ for the decay H ^ ZZ. To one- loop accuracy A^[pR can be 



'-W'^W i=0,l 



„2 2 -"-'00 _ 



yn^^^^ _ - -is- Sp)J:'^{ 

S Sp 



(35) 



Owing to Sp = — iMnrn + C(q;^), the right-hand side of (|35D differs from (p6[) and 
(p8|) only by higher-order contributions. 

Note that the freedom in splitting the matrix element into SMEs and form factors 
is also present in this approach, i.e. in (^) we could also introduce functions fij{s), 
as described in the first part of this section. In the considered case one could avoid 
this ambiguity by expanding also the SMEs occuring in Ai^'^'^ about Sp. However, in 
more complicated situations it is not always possible to include the wave functions in 
the pole expansion. If one then excludes the SMEs from the pole expansion, as for 



instance advocated in Ref. [ one is again confronted with the problem of violating 
Ward identities. 

5.2 Dyson summation within the background-field method 

A different approach to introduce a finite width near resonances is to Dyson-sum the 
self-energy corrections. It is a well-known fact that in the conventional formalism Ward 
identities, which in particular rule the gauge cancellations, are violated if Dyson summa- 
tion is applied. However, in Ref. [ |13[ it has been shown that Dyson summation within the 
background-field method (BFM) (see Ref. [ ^ and references therein) does not violate 
the Ward identities if all one-particle-irreducible corrections are taken into account in the 
same loop order. Dyson summation naturally arranges the reducible parts of amplitudes 
in a way that results from forming trees with vertex functions joined by full propagators 
(inverse two-point functions). For the process under consideration this simply amounts 
to writing the (one-particle-reducible) lowest-order, self-energy, and vertex contributions 
in the following way: 

A^Born + ^A^self + (^Al vertex 

- -MiPR - ^ ^ M,^ r-Ml + T.^{r) ' ^^^^ 



CwS 



W'W i,j=0,l r=s,t,u 



where T,^ and Ff^^ denote the renormalized self-energy and form factors, respectively, 
in the BFM and 

-^S? = {erks){e2-h){el-h){el-h)/s^ M'tf = {erk,){e2-h){el-k2){el-h) / s\ (37) 

The s-channel part of ( PB[ ) is formally identical to (PTf). Note, however, that we use 
(^) in the following with form factors and self-energy in finite, i.e. one-loop, order of 
perturbation theory. 
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The complete one-loop matrix element is obtained by adding the (one-particle-irredu- 
cible) box contributions and multiplying everything with the (UV-finite) wave-function 
renormalization factor \fRz for each external Z boson, 

A^BFM = (A^IPR + A^box) [Rzf . (38) 

Since the wave-function renormalization factor Rz multiplies the complete matrix element, 
we can simply use its 0{q) approximation 

{Rzf = 1-2 ReS'|^(M|). (39) 

The matrix element ( pSj ) fulfills all relevant Ward identities, and, in particular, does 
not violate unitarity at high energies. However, it depends on a gauge parameter via 
higher-order (at least two-loop) corrections which are not completely taken into account. 
This is nothing but a result of the fact that a Dyson summation is always arbitrary to some 
extent. Note also that the on-resonance self-energy is unique and equal to the physical 
quantity iMnTH. 

In the conventional formalism the matrix element after Dyson summation depends 
not only on the gauge but in addition on the choice of the field renormalization. In 
the BFM the matrix element is actually independent of the field renormalization. This 
can be seen as follows: The field renormalization is fixed by background field gauge 
invariance up to a UV-finite linear transformation of the renormalized fields. Such a 
linear transformation turns the linear Ward identities for the background-field vertex 
functions into Ward idenities for transformed vertex functions with the same structure. 
These modified Ward identities are still exactly valid even for full, i.e. Dyson-summed, 
propagators. However, the effects of the linear transformation cancel in ^'-matrix elements, 
thus giving a unique answer.' 

The resonant part of the Dyson-summed one-loop matrix element (|36|) in the BFM 
reads for s = My 



M 



BFM 



This differs from the resonant part of the pole-scheme amplitude (|30|), 

^, e^Ml - S'^(M^)) + E^=oAM\^ + M^^hFF"{Ml) 



M 



pole 



A/2 iMnrn 



(41) 



in 0{a) in relative terms, in accordance with the discussion after (|27|). Moreover, (|40|) 
is even gauge-dependent in this order on resonance, whereas ( ^T] ) is manifestly gauge- 
independent. The bulk of these effects can be attributed to the contribution of S''^(M^) 
and thus to the different wave-function renormalization in the BFM. 

In order to obtain the cross-section on resonance also in 0{a) accuracy the imaginary 
part of the Higgs-boson self-energy has to be included in two-loop order. In the pole 
scheme this is equivalent to the introduction of the (9(Q;)-corrected Higgs-boson width in 
the propagator. However, in the BFM approach all two-loop corrections are required in 
order to preserve the Ward identities. 
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Figure 10: Lowest-order diagrams for XX 



XX 



6 Cross-section for longitudinal Z bosons from the equivalence theorem 



The corrections of O (aM^/s^M^) to longitudinal gauge-boson scattering processes 
have been calculated in the literature [ H, ^ || using the equivalence theorem (ET) 
[ ^, [TD|, [TT|, |r^. The ET relates amplitudes involving longitudinal gauge bosons 
with those involving the associated would-be Goldstone bosons in the high-energy limit. 
Because the latter amplitudes are much easier to be calculated, the ET was frequently 
used to obtain cross-sections in the high-energy limit. 



6.1 Equivalence theorem within the background- field method 

When using the ET in higher-order calculations one has to be careful to include all 



correction factors that result from renormalization and amputation [ 11|. It has been found 
in Ref. [ Q that this is particularly easy within the BFM. In this formalism the matrix 



elements for external longitudinal vector bosons are directly obtained from the amputated 
Green functions with the corresponding would-be Goldstone boson fields multiplied with 
the wave-function renormalization constants of the gauge bosons. Moreover, in contrast 
to the conventional formalism, in the BFM the ET is valid with and without Dyson 
summation. 

6.2 Results for ZlZl ZlZl 

We want to apply the ET to the process ZlZl ZlZl in the framework of the BFM 
and investigate the accuracy of the corresponding predictions. 

To this end we have to consider the process xX XXi where x is the would-be 
Goldstone boson associated with the Z boson. In lowest order the four diagrams of 
Fig. ITOl yield 



A A XX *XX 



'Ml 



'Born 



As 



w'-w 



3 + 



Ml 



Ml 



+ 



Ml 



Ml 



+ 



Ml 



u 



Ml 



(42) 



For s ^ M| the matrix element for ZlZl — > ZlZl approaches the one for xX ~^ XX fo^' 



any value of Mh [ 12|. We note that this agreement is destroyed if one includes the finite 
Higgs-boson width in the way discussed in Section ^3] without expanding also the SMEs 
or appropriately adjusting the function /(s) in (p^). 

Following the treatment of Section |5.2| , the matrix element for xX ~^ XX including 
0{a) corrections is given similarly to (^) by 



M 



xx-^xx 



A A XX ^XX 
J^l 1 



4PR 



4^2 



Ml 



A A XX *XX 



^box 



(43) 
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with 

Because A^^p^'^'*^ vanishes for s 3> as M^/ s, the boxes dominate the matrix element 
for XX ~^ XX ci't high energies. The matrix elements for XX ~^ XX ^-re calculated in ex- 
actly the same way as for ZZ — >■ ZZ. In contrast to ZZ TIL the one-loop corrections 
to XX ~^ XX possess a transparent form; they are explicitly presented in Appendix p. 
We recall that the usage of the BFM is crucial for the ET to work for Dyson-summed 
amplitudes. Because the would-be Goldstone bosons are scalars, no polarization vectors 
and no SMEs occur, and no unitarity cancellations between individual contributions take 
place. This simplifies the calculation considerably. 

The Higgs-boson resonance can be treated exactly as for ZZ TIL. Since the wave 
functions are trivial constants, and thus no split into SMEs and invariant functions is 
necessary, the ambiguity in applying the pole expansion is absent. 

In order to improve the accuracy of a calculation via the ET one can combine the full 
lowest-order matrix element with the 0[oi) corrections from xX ~^ XX resulting in 

Admixed = A^L'r'^'^ + ^MW^""- (45) 

This treatment is, however, not possible if one uses Dyson summation, because in this 
case the lowest-order matrix element cannot be linearly separated from the one-loop 
corrections. 

6.3 Heavy-Higgs-boson effects 

In the literature an approximation for the matrix element TW-^-^^^^ by the leading 
contributions for s, M\ ^ M| was frequently used [ |^, |, ||, ||, ^, §] . In this approximation 
the building blocks of p3| ) and (0) take a particularly simple form. In the BFM we find 
in this limit 



" ^|:^{3M|[35o(r, Mh, Mh) + 5o(r, 0, 0) 



47r 8s2^Mw 



3fio(M2 , Mh, Mh) - Re{5o(M2 , 0, 0)}]} + bZ^{r - Ml 



pxxB^^^ = --^-^%^{2M2 [3Co(r, 0, 0, Mh, Mh, 0) + Co(r, 0, 0, 0, 0, Mh)] 
+ 3So(r, Mh, Mh) + 5So(r, 0, 0) + 45o(0, 0, Mh) 

- 9i?o(M2 , Mh, Mh) - 3 Re{5o(M2, 0, 0)} - ^} + Uz^ + bZ^, 

Ml™ = a'—^\-^{2Ml[Do{0, 0, 0, 0, s, t, Mh, 0, Mh, 0) 

+ Do(0,0,0,0,s,t,0,MH,0,MH)] 

+ [Co(s, 0, 0, Mh, Mh, 0) + 3Co(s, 0, 0, 0, 0, Mh)] 

+ 5o(s,Mh,Mh) + 115o(s,0,0) 

- 9Bo{Ml, Mh, Mh) - 3 Re{Bo{Ml, 0, 0)} - 1 
+ {s-^t,t-^u) + {s u,t ^ s)} + 26Z^, 
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Rz = 1, 



(46) 



where Bq, Cq, and Dq are scalar one-loop functions 
malization constants read in the BFM 



[ |18|, . The wave-function renor- 



5Zfj = 5Z^ 



a 



Ml 



47r8s^M2,- 



(47) 



The above results are in agreement with those of Ref . [ |^ . Since the approximation 
(^) merely involves corrections of (9 (aM^/s^M-^), it follows from power counting 
(see Ref. [ 0] and references therein) that only diagrams with internal scalar lines con- 
tribute. We note in passing that the terms of O (aM^/s^M^) originate entirely from the 
SU(2) sector of the SM, i.e. they could also be obtained from the corresponding reaction 
W^W^ W^W^ in the pure SU(2) gauge theory. 

We have checked that in the limit M| ^ s ^ Ml the amplitude reduces to 



E 



r=s,t,u ISS^M^ 



In 



Ml 



3^3: 



TT 



-r — le 
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y 



(45 



as already given in Refs. In this context we remark that the result (P8| ) can be 

most easily obtained from the general structure of the heavy-Higgs-boson limit of the SM. 
The matrix element JV[^^^^^ ( ^8] ) gets contributions only from £4 and £5 of the effective 
Lagrangian of Ref. [0], which quantifies the heavy Higgs-boson effects, and from the 
(three) irreducible graphs in the gauged non-linear cr-model which contain only quartic 
scalar couplings. 



7 Discussion of numerical results 



7.1 Computational details 



For the calculations we use the following parameter set [ ^ 



a-^ = 137.0359895, 
= 0.51099906 MeV. 



= 105.658389 MeV; 



rrir 



1771.1 MeV. 



M7. 



rrir 



-- 91.188 GeV, 

47.0 MeV , 
1.55 GeV, 
180 GeV , 



Mw = 80.26 GeV , 

47.0 MeV, 
150 MeV , 
4.5 GeV . 



ma 
rris 



(49) 



The masses of the light quarks are adjusted such that the experimentally measured 
hadronic vacuum polarization is reproduced [ For a Higgs-boson with a mass of 
Mh = 700 GeV these parameters yield the lowest-order decay width Fh ~ 175.29 GeV, 
i.e. about one fourth of the mass. 

The various independent calculations described in Section ^?T] agree numerically typ- 
ically to ~ 10 digits apart from the regions close to the boundaries of phase space. At 
these boundaries the reduction of tensor integrals to scalar integrals breaks down. We 
avoid these regions by using the angular cut ^cut = 10°, which also removes the Landau 
singularities in the fermionic boxes for energies above about 500 GeV. 
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7.2 Corrected cross-sections 



The integrated cross-sections for unpolarized, purely transverse, and purely longitudi- 
nal Z bosons in lowest order and including the one-loop corrections are shown in Figs. ITT 



and |T2| for Mh = 100 GeV and Mh = 700 GeV, respectively (repeating information from 
Figs. @, ^ ^ and H). In the case of Mh = 100 GeV no finite Higgs-boson width is in- 
troduced; for Mh = 700 GeV we apply Dyson summation within the BFM using the 
renormalization scheme of Ref. [|16|. The Higgs-boson- mass dependence of the cross- 
section for purely transverse Z bosons is below 10% including the Higgs-boson-resonance 
effects. The corresponding lowest-order cross-section is very small at high energies and 
not visible in Fig. 



The differential cross-sections for various energies are shown in Figs. |I3] and |1J. For 
small energies the angular variation of the cross-sections stays within one order of mag- 
nitude. For high energies and a small Higgs-boson mass the corrected cross-sections are 
strongly peaked in the forward and backward directions, while the lowest-order cross- 
sections are relatively flat in the considered angular region. The cross-section for purely 



longitudinal gauge bosons has kinematical zeros 

if Mh < V 1 + V^Mz ^ 150 GeV, which 
move towards the forward and backward directions with increasing energy. For a large 
Higgs-boson mass the cross-section for purely longitudinal Z bosons, which dominates in 
this regime, becomes flat and therefore also the unpolarized cross-section. 

7.3 Higgs-boson resonance 



In Figs. |T^ and |T6| we compare several different treatments of the Higgs-boson reso- 
nance using Mh = 700 GeV. We include the lowest-order (Born, Fh = 0) and the corrected 
(full, Fh = 0) cross-sections for vanishing Higgs-boson width for reference. We show the 
lowest-order (Born, pole scheme) and the corrected (full, pole scheme) cross-section in the 
pole-scheme treatment given in (EBf) and (EBf) as well as the corrected cross-section in the 



modified pole scheme (full, mod. pole scheme) according to (|30|) with f{s) = M^/s. In 



addition we give the cross-sections resulting from Dyson summation according to (130) and 
(p8|) in the BFM (BFM, Dyson) and the corresponding one in the conventional formalism 
(conv., Dyson). Apart from the Dyson-summed cross- sect ions all others are identical in 
the BFM and in the conventional formalism. Since the unpolarized cross-section is dom- 



inated by the one for purely longitudinal Z bosons for a large Higgs-boson mass. Fig. |15 
holds essentially also for the latter cross-section after multiplying by a factor nine. 



The crucial differences between the various treatments can already be seen in Fig. 
which shows the integrated cross-section. Owing to the crude resolution the pole-scheme 
cross-sections with or without 0{a) corrections cannot be separated from the correspond- 
ing cross-sections for Fh = at high energies. The deviation of the Dyson-summed BFM 
cross-section is due to higher-order corrections that become increasingly important with 
energy. The Dyson-summed conventional cross-section deviates more for energies above a 
few TeV and becomes completely wrong for energies higher than 10 TeV. This results from 
the violation of the Ward identities which leads to unitarity violation at high energies. 

The differences between the various treatments of the Higgs-boson resonance can be 
seen more clearly in Fig. 0, where the corrected cross-sections are shown normalized to 
the one in the modified pole scheme. The difference between the pole scheme and the 
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Figure 11: Lowest-order and corrected integrated cross-sections for various polarizations 
at Mh = 100 GeV 
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Figure 12: Lowest-order and corrected integrated cross-sections for various polarizations 
at Mh = 700 GeV (The lowest-order TTTT cross-section is not visible.) 
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Figure 13: Lowest-order and corrected differential cross-sections for various polarizations 
and CMS energies at Mh = 100 GeV 
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Figure 14: Lowest-order and corrected differential cross-sections for various polarizations 
(as indicated in Fig. |13D and CMS energies at Mh = 700 GeV 
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Figure 15: Integrated unpolarized cross-section at Mh = 700 GeV for various treatments 
of the Higgs-boson resonance 
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Figure 16: Relative deviation of various treatments of the Higgs-boson resonance from the 
modified pole-scheme result for the integrated unpolarized cross-section at Mh = 700 GeV 
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modified pole scheme is below 2% and becomes small at high energies. Note, however, 
that by using dimensionful SMEs the pole-scheme cross-sections could become completely 
wrong at high energies owing to spurious unitarity-violating terms. 

In the resonance region the Dyson-summed cross-sections deviate from the cross- 
sections in the modified pole scheme by up to 19% and 7% in the BFM and the con- 
ventional formalism, respectively. This difference is due to the fact that our calculation 
near the resonance is only of 0(1) accuracy since the lowest-order contribution in the res- 



onance denominator vanishes on resonance (cf. Section ^.1|) . The size of these differences 
and the correction of 24% of the pole scheme calculation on resonance set the typical scale 
for the missing 0{a) corrections in the resonance region. 

7.4 The cross-section for ZlZl — > ZlZl and the equivalence theorem 

Finally, we want to investigate the numerical accuracy of the ET. We distinguish the 



cases without and with a Higgs-boson resonance. In Fig. |1^ we consider the case of no 
Higgs-boson resonance (Mh = 100 GeV). We show the lowest-order cross-section (Born) 
calculated from the ET normahzed to the lowest-order cross-section for ZlZl —>■ ZlZl, the 
fully corrected cross-section (full) calculated from the ET and the cross-section obtained 
from the matrix element ( ^3]) (mixed) both normalized to the fully corrected cross-section 
for ZlZl ^ ZlZl. The quality of the ET at Ecus = 1 TeV (2TeV) is about 17% (5%) 
for the lowest order, 24% (10%) for one-loop, and 6% (4%) for one-loop mixed. As 
expected, the one-loop mixed approximation is substantially better than the simple ET 
cross-section. 

In Fig. |l^ we investigate the accuracy of the ET in the presence of a Higgs-boson 
resonance at Mh = 700 GeV. We show again the ratios of the lowest-order and corrected 
cross-sections obtained using the ET and from the direct calculation. In the lowest- 
order cross-section we include the finite width naively (Born) and in the modified pole 
scheme (mod. Born) ([29|). The lowest-order cross-section from the ET approaches that 
of the modified pole scheme at high energies. Including the finite width naively leads to 
a cross-section that deviates at high energies from these two cross-sections by a factor 
l + T^/9Ml ^ 1.007 for Mh = 700 GeV (cf. Section For the corrected cross-sections 



we have applied Dyson summation within the BFM. Because of the Dyson summation 
the mixed case does not make sense anymore. Instead we show the cross-section resulting 
from the O (aMn/sw-^w) approximation of the RCs (^61) normalized to the fully corrected 
cross-section for ZlZl — * ZlZl. The ET works much better for a heavy Higgs boson. At 
Eqms = 1 TeV (2 TeV) we now find a deviation of 8% (2%) for the lowest order [using 
the modified Born according to (^^] and 6% (2%) for the corrected cross-section. For 
energies above 2 TeV the deviation between the corrected cross-sections is practically 
equal to the deviation between the lowest-order cross-sections. The O {aM^/ s'^M^) 
approximation ( ^61) works well in the regime M| -C M^ <^ s <^ Mh/M|, where the upper 
limit for the energy results from the neglect of corrections proportional to s/M| with 
respect to the ones proportional to M^/M|. For Mh = 700 GeV this restricts the energy 
to ^/s ~ 1 — 3 TeV, which is nicely reflected in the figure. 
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Figure 17: Relative deviations of the ET predictions for ZlZl 
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8 Conclusions 



Owing to the strong sensitivity to the gauge-boson and scalar self-interactions, scat- 
tering of massive gauge bosons found continuous interest in the literature, where the em- 
phasis was directed to strong-coupling effects for longitudinally polarized gauge bosons. 
We have supplemented the existing results for the enhanced radiative corrections of order 
O («M^/s^M^) by the complete 0{a) corrections to ZZ TIL for arbitrarily polarized 
Z bosons. 

At high energies the radiative corrections are found to be large, at several TeV they 
are typically of the order of the lowest-order cross-sections. Whereas the cross-section 
for purely transverse Z bosons at high energies is totally negligible in lowest order, the 
corrections enhance this cross-section such that it becomes one of the dominating channels. 

The introduction of a finite Higgs-boson width in order to describe the resonance well 
is a non-trivial task. We have compared different approaches, viz. different variants of 
the Laurent expansion about the complex pole and the Dyson summation of self-energies, 
where the latter has been performed both in the conventional formalism as well as in 
the background-field formalism. From a theoretical point of view, the background-field 
approach is the most convincing one, since it naturally guarantees a reasonable cross- 
section also far above the resonance, where the vahdity of Ward identities is crucial to 
imply the necessary gauge cancellations. However, in order to obtain a relative precision 
of Oici) on resonance one would have to perform a complete two- loop calculation. In 
order to obtain the same precision on resonance in the pole scheme only the imaginary 
part of the self-energy has to be evaluated at two loops. However, as the pole scheme 
and the other mentioned methods do not care about the Ward identities, theoretical 
uncertainties may get out of control in the presence of gauge cancellations. Using the 
pole scheme carelessly can lead to unitarity-violating terms at high energies, and Dyson 
summation within the conventional formalism in fact yields a totally wrong cross-section 
in the high-energy limit. 

We have investigated longitudinal Z-boson scattering ZlZl — > ZlZl in more detail and 
performed a complete Oipi) calculation using the Goldstone-boson equivalence theorem. 
For a center-of-mass energy of 1 TcV (2 TcV) the deviation of the equivalence theorem 
from the exact 0{ol) result is about 24% (10%) and 6% (2%) for a Higgs-boson mass Mh of 
100 GeV and 700 GeV, respectively, with an asymptotic approach in the high-energy limit. 
The frequently used approximation by the enhanced corrections of O (aM^/s^M-^) for 
a heavy Higgs boson is good for energies of a few TeV but gets worse with increasing 
energy. 

Although Z-boson scattering is the simplest representative of massive gauge-boson 
scattering, it contains the typical features such as the Higgs-boson resonance and enhanced 
heavy-Higgs-boson corrections. In contrast to other gauge-boson scattering processes, 
the lowest-order cross-sections for transverse Z-boson scattering are suppressed and no 
real photon radiation needs to be considered in ZZ — > ZZ. Nevertheless, we expect that 
our results at least qualitatively carry over to the other massive gauge-boson scattering 
reactions. 
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Figure 19: Generic form of some singular contributions to ZZ 'Ziff, which are related 
by cutting rules. 

Appendix 

A Discussion of the Landau singularity in box diagrams 



In Section |4.4| we have briefly discussed the Landau singularity which occurs in some 
fermionic box diagrams if the fermion mass m fulfills m < Mz/2. Although we have argued 
that this singularity is unphysical and only caused by the use of the equivalent vector- 
boson approximation for pp or ee collisions, it is nevertheless interesting to investigate 
some formal properties of the singularity. 

From general considerations (e.g. about unitarity) one expects that the singularity 
drops out in the fully inclusive cross-section, i.e. if all possible final states are taken 
into account. We have verified this compensation by explicitly calculating the singular 



contributions of Fig. 19 to the inclusive cross-section ZZ 4/. The shaded circle in 



Fig. |I9| represents any regular graph for ZZ ZZ, i.e. only the cuts that are explicitly 



shown in Fig. |T9| are relevant for the singular contributions. If the shaded circle also 



contains the singularity more cuts have to be considered. We restrict ourselves to the 

V3 = /c| = M|. 



case where the produced Z bosons are on their mass shell, k\ = k\ = M| 



The singularity in the loop integral of the upper graph of Fig. |1^ stems from integration 
momenta gi ~ g with 

= -^{ki - hr + ih + k^r, = — 1^ — , (Ai) 
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and occurs for (f m^, which follows from the Landau equations for the integral. Thus, 
the singular contribution of the virtual graph is simply obtained by setting gi — g in 
the numerator such that the remaining integral is proportional to the Dq function, the 
singular contribution of which is given in (|2T|) . 



The singular contributions in the lower graphs of Fig. occur in the phase-space 
integrals of the produced fermion-antifermion pair and can be obtained from the corre- 
sponding scalar integrals analogously to the loop integral. For instance the relevant scalar 
integral for the lower left graph is given by 

'^'^ ' {2'nY5^'Kh-k-k')- (A2) 



(27r)32A;0 J {2n)^2k'0 ' ' ' " ' {qf - + ie)(g| - + ie 

where 

qi = k — ki, q4 = k — ki + ki = k2 — k', k"^ = k'"^ = vr? . (A3) 

Inspection of the Landau equations for / reveals that the singularity in / originates from a 
point in phase space where the fermion momenta k and k' are coplanar with the scattering 
plane spanned by the k;. Explicit calculation yields 

rising = ~Y6^ ' (^^) 



which is also valid for the corresponding scalar integral for the lower right graph of Fig. |l£ . 

The relation between the singular contributions of and / guarantees the cancellation 
of the singularity in the sum of the graphs of Fig. |19| if the produced Z bosons are on 
shell. However, the cancellation in general is incomplete if individual flavors or spins in 
the fermionic final state are observed or if phase-space cuts are applied. 

Therefore, we conclude that a careful analysis of the actual physical realization of the 
underlying process is mandatory if such singularities appear for physical situations. For 
the subprocess ZZ — >■ ZZ this means that one has to go one step back and to consider 
the full reaction including the production mechanism of the incoming Z bosons in more 
detail. 



B One-loop corrections to XX ~^ XX 



In Section 3.3 we have explicitly given the leading corrections to longitudinal Z-boson 



scattering in the limit s, ^ M|. More generally, the complete corrections to xX 
also take a relatively simple form in contrast to the formulae for ZlZl ZlZl, 



XX 



a 



47r 8s^Mw 



{9M^5o(r,MH,MH) 



'w 



+ 2(M^ + m'^M^ + 12M^ - 8rM^)5o(r, Mw, M. 
+ (M^ + 4M^M| + 12M| - 8rMi)Bo{r, Mz, Mz) 
+ 3M^Eo(0, 0, Mh) + 2M^(M2 + 6M^)5o(0, 0, Mw) 
+ M|(M^ + 6M|)5o(0, 0, Mz) - 24M^ - 12M| 

/ 



(r - Amj)B, 



o{r,mf,mf) -2mjBo{0,0,mf) | 
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pxx 



H/ 



5Ml + 5Zfj{T-Ml), 
1 



a 



{ 



6M 



H 



(M^ - M|)^ + 2M|(r - 2M|)J Co(r, M|, M|, Mh, Mh, Mz) 
+ 2(M^ - 3M^Ml - 6M| + ArM^M^)Co{r, M|, M|, Mz, Mz, Mh) 
+ 8M2,(M2 + 4M2,)(r - 2M|)Co(r, M|, M|, Mw, Mw, Mw) 
+ 3M^(M^ + 2M|)So(r, Mh, Mh) 
+ 2(M^ + 4M2M^ + 12M^)So(r, Mw, Mw) 
+ (3M^ + GM^Ml + 8Ml)Bo{r, Mz, Mz) 
+ 4(M^ - MifBoiMl Mh, Mz) - 16M^ - 8M^ 
+ SiV^m} [(2M| - r)Co(r, M|, M|, m/, m/, m/) - 2Bo{r, ruf, m/)] } 



9„2 



2M2, 



Mi 



5t 



2sw MwMi 



'box 



= a 



{2 [(Mi 



X Do(M|, M|, Ml Ml s, t, Mh, Mz, Mh, Mz) 



+ 2 (Mi 



M|)2 + 2M|(t-2M|) 



X Do{Ml Ml Ml Ml t, s, Mh, Mz, Mh, Mz) 
+ 16M^ [{s - 2Mlf + (t- 2Mlf 

X DoiMl Ml Ml Ml s, t, Mw, Mw, Mw, Mw) 
+ 4(Mi + 2M|) [(Mi - Mif + 2M|(s - 2M|)" 

xa3(s,M|,A^,MH,MH,Mz) 
+ 4 [SM^iM^ - Mlf + 4M^(s - 2M|)] Co{s, Ml Ml Mz, Mz, Mh) 

+ 16M^(Mi + 4M^)(s - 2M|)Co(s, M|, M|, Mw, Mw, Mw) 

+ (Mi + 2M|)25o(5,MH,MH) 

+ 2(Mi + 4MiM2, + 12M^)So(s, Mw, Mw) 

+ (9Mi + 8M|)5o(s, Mz, Mz) - IGM^, - 8M^ 

+ Y.^N]^) [{st - 2Ml)Do{Ml Ml Ml Ml s, t,mf,mf,mf,mf) 
f 

+ 4(2M| - s)Co(s, Ml Mlmf,mf, ruf) - 4So(s, m/, ruf) 

+ (s ^ t,t ^ u) + (s — >• M, t — >• s) I 

e 5t 



+ 2SZ,-^-^ 



SM^^SM^^ 



Mi 



2sw MwMi 



+ 2SZ. 



X' 



(Bl) 



32 



where the sum over / extends over all fermion flavours, and Nj denotes the color factor 
for the fermion /. The scalar four-point function is defined as {pi + P2 + + P4 = 0) 



Doipl,pl,pl,pl, (Pi +^'2)^ ip2+P3y, mi, 1712, 1713,1714) 
d\ 1 



(B2) 



ivr^ (g^ - ml) [{q + pi)^ - mj] [{q + Pi + ^2)^ - "^2] [il - PiY - "^2] 



Note that these results are derived within the BFM and include the heavy-Higgs-boson 
corrections of ( ^6]) as special case. For the sake of simplicity the explicit expressions for 
the counterterms are left open; they are easily calculated in the renormalization scheme 
of Ref . \ M 
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